arXiv:hep-lat/9706020v2 15 Nov 1997 


On the SU(2)-Higgs Phase Transition 


Isabel Campos * 
February 7, 2008 


Departamento de Fisica Teorica, Facultad de Ciencias, 
Universidad de Zaragoza, 50009 Zaragoza, Spain 
e-mail: isabel@sol.uiiizar.es 


Abstract 

The properties of the Confinement-Higgs phase transition in the 
SU(2)-Higgs model with fixed modulus are investigated. We show that 
the system exhibits a transient behavior up to L=24 along which, the 
order of the phase transition cannot be discerned. To get stronger 
conclusions about this point, without going to prohibitive large lattice 
sizes, we have introduced a second (next-to-nearest neighbors) gauge- 
Higgs coupling (ac 2 ). On this extended parameter space we find a 
line of phase transitions which become increasely weaker as K 2 —> 0. 
The results point to a first order character for the transition with the 
standard action (k 2 = 0). 


* Supported by a Fellowship from Ministerio de Educacion y Cultura, Spain 


1 



1 Introduction 


The generation of mass in the electroweak sector of the Standard Model 
(SM) |]l| relies on the Higgs mechanism. Because of this fact large amount 
of work has been spent, both, perturbative and non-perturbatively in order 
to understand the continuum limit of gauge-Higgs models. Underlying is the 
very question: Does a non-trivial QFT exists in the limit of inhnite cut-off ? 
(see [@] for a review). 

On the one hand, it is almost rigorously proved p, Q that the pure scalar 
sector, namely the A<h^ model, leads to a trivial theory when the cut-off 
is removed. Of course this is an academic model, and the question is if 
the coupled gauge-scalar model could produce a non-trivial theory in the 
continuum. 

Perturbatively, the answer seems to be negative, only theories asymptot¬ 
ically free in all the couplings can be constructed But realistic theories 
have at least one coupling which is not asymptotically free. A quite accepted 
scenario describes the SM as an effective theory, with a hnite cut-off above 
which the theory is no longer valid. Within this approach an upper-bound 
for the Higgs mass can be calculated (see 0 for a review). 

However, the non-perturbative sector of the SM is not yet completely 
understood. The strong self-coupling allowed for the Higgs held, renders 
perturbation theory useless, and one is forced to use non-perturbative meth¬ 
ods to get insight into the properties of the model in this region of couplings. 

The Higgs sector of the SM can be approximated by the SU(2)®U(1) 
Higgs model. This approximation is expected to behave reasonably well for 
the Yukawa couplings between the fermions and the Higgs held, excepting the 
coupling of the top quark, are small. Also, since the U(l) and the SU(2) gauge 
couplings are related through the Weinberg angle by ~ 0.27(7su(2), we 
can start with the SU(2)-Higgs model, neglecting the U(l) degree of freedom, 
as a hrst approximation to describe the electroweak interaction. 

In particular, we shall be interested in the limit in which the modulus 
of the Higgs held is frozen. In the usual notation it correspond to the limit 
A = cxo, being A the parameter controlling the radial degree of freedom of the 
Higgs held. The phase diagram for this model is well known 0JJ- A phase 
transition (PT) line separates a region where the scalar particles are conhned 
in bounded states (conhned phase), from another region where the symmetry 
SU(2)®SU(2) is spontaneously broken, and the spectrum consists of the W’s 
gauge bosons and the Higgs particle (Higgs phase). In this sense SU(2)-Higgs 
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is similar to QCD: there is a phase in which gauge color helds form glueballs 
and quarks are conhned into hadrons, and another phase characterized by 
the onset of the Higgs mechanism. 

The PT line ends at some hnite point of the parameter space being both 
regions analytically connected. So, strictly speaking, we should talk about a 
single phase, the Conhned-Higgs phase (see Figure]^). 
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Figure 1: Schematic representation of the phase diagram of the SU(2)-Higgs model 
with frozen modulus (A = oo). 


In the scaling region, this model has been extensively studied for small 


tinctly hrst order. The PT weakens when increasing A. In the limit A = cxo it 
is generally believed that the transition is still hrst order, though extremely 
weak. However, in our opinion, a higher statistic study in the limit A = cx) 
is still lacking since the results up to date are not conclusive, as we shall 
demonstrate. To shed some light on this problem we have studied the model 
in an extended parameter space too. However, our motivation is not only 
performing such a large statistics study, but also we want to extract gen¬ 
eral properties of weak hrst order phase transitions in coupled gauge-Higgs 
systems. With this purpose, we have added an extra positive gauge-Higgs 
coupling between next-to-nearest neighbors to the standard action. We shall 


m and intermediate 0,0, m values of A, where the transition is dis¬ 
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use it as a parameter to study the weakening of the PT when this extra 
coupling is tuned to zero. 

In the next section we describe the model, and summarize previous re¬ 
sults. In section 3 the numerical method and analysis techniques are detailed. 
Section 4 contains the results. Finally, the last section is devoted to conclu¬ 
sions. 


2 The Model 

The SU(2) lattice gauge model coupled to an scalar held, in the fundamental 
representation of the gauge held can be described by the action 


S, = /3^|l - iTl'f/J - h, + ^L) + 


X,(l 


+X ^[<h^(a;)<F(x) — 1]^ -I- ^ <h^(a;)<F(a;) 


( 1 ) 


Where U^(a;) represents the link variables, and Up are their products 
along all the positive oriented plaquettes of a four-dimensional lattice of side 
L. 

The scalar held at the site x is denoted by <F(a;), being A the parameter 
controlling its radial mode. In the limit (A = oo, (3 = cxo) the model becomes 
a pure 0(4)-symmetric scalar model. 

As we pointed out in the previous section, the PT line ends at some hnite 
value of the parameters (/3, ki) The endpoint moves towards larger /3 values 
as A increases. For A > 0.1 the endpoint crosses to the /3 > 0 region, and 
in the limit X = oo the phase transition ends at (ki ~ 0.6, (3 ~ 1.6) [^. It 
is commonly believed that the transition at this point is second order with 
classical critical exponents [II^, however a careful numerical study would be 


necessary. 

In the scaling region, and for hnite A, the phase transition turns out to 
be hrst order. Also, the transition becomes weaker as /? or A increases. In 
particular in the limit j3 = oo (spin model) the transition is second order 
with classical critical exponents [p!4[ . 

In the limit A = cx) the situation is less transparent, and it is not clear 
whether the transition is weak hrst order or higher order. 
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The study of the model with the parameter A = oo is equivalent to fixing 
the modulus of the Higgs held, $1 ■ $ = 1. In the pioneer work, ||^, the PT 
was considered hrst order for hnite A and second order for A = cx). Later, 
larger statistics, and the hypothesis of a universal behavior of the PT for all 
values of A, seem to point to a hrst order character of the PT in the scaling 
region. Nowadays, though it is generally believed that the transition is still 
hrst order in this limit, the numerical proofs BUI on which rely these 
statements are not conclusive, as the authors safely conclude, because the 
statistic and lattice sizes are not enough for excluding the possibility of a 
higher order phase transition in the edge A = cx). 

The model with hxed modulus is described by the action 

^oo = - ^Trt/p] - + (2) 


As we shall show below, we have simulated this model up to lattices T=24, 
and the result is still compatible with a second or higher order PT, however, 
we have no indications of asymptoticity in the behavior of the observables, 
and the results are compatible with a very weak hrst order PT too. This 
means that larger lattices are needed to overcome these transient effects, but 
the added difficulty here is that such lattices would suffer of termalization 
problems, and severe autocorrelation times. Altogether makes this approach 
too CPU expensive for nowadays computers. 

In order to get a more conclusive answer without going to prohibitive large 
lattice sizes, we have studied the model in an extended parameter space. For 
this purpose we have introduced a second coupling between the gauge and 
the scalar held connecting next-to-nearest neighbors on the lattice, in such a 
way that the new action reads 

-S'= S'oo-^K2 ^ Tr^\x)[U^{x)Uy{x+^ji)+Uy{x)U^{x+v)]^{x+^+v) (3) 


Within this parameter space we expect to get a global vision on what the 
properties of the PT are, and also, to give a stronger conclusion about the 
order. In the region of K 2 positive, (competing interaction ehects appear if 
K 2 < 0) this extended model is expected to belong to the same universal¬ 
ity class than the standard one (^2 = 0) since both models posses the same 
symmetries. The ehect of the new coupling ^2 > 0 is to reinforce the transi¬ 
tion but should not change the order if the system has not tricritical points. 
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An example of such behavior appears in the 0(4)-symmetric a model with 
second neighbors coupling. This model presents a ^ 2 ) line of phase tran¬ 
sitions which is second order, since the model with ^2 = 0 shows a second 
order PT too [pT[| . 

The action (^) has the following symmetries: 

• = ^2 = 0 . 


(3 ^ -f3 
U^{x) ^ 


X) 


> — Ki, <h(x) hxed 


U^{x) ^ -U^{x) 


—i> —Hi, Ufj,{x) hxed 


<h(x) ^ 


(4) 

(5) 

( 6 ) 


The action is not symmetric under the change K 2 —>■ — ^ 2 - The existence 
of couplings Hi and H 2 with opposite signs would make frustration to appear, 
and very different vacua are possible |]^. This region present problems 


when one tries to implement rehection positivity, however, the possibility of 
dehning a continuum limit at this region is not discarded a priori [|T^. In 


this work we are interested in the regions free of frustration effects. Taking 
into account the symmetry properties of the action, the phase diagram in the 
region K 2 > 0 will be symmetric with respect to the axis Hi = 0, and hence, 
we can restrict the study to the quadrant Hi > 0. 

We dehne the normalized energy associated to the plaquette term 


Eo = — V(1 - -TtU„) 

Wo ^ 2 

'J p 

and also the energies associated to the links 

Ei = ^Yl ^Tr$^(a:)f/^(a:)<h(a: + /i) 


(7) 


( 8 ) 


W = ^ V ]TT^Hx)[Up{x)U,{x + ^i) + U,.{x)U^{x + ly)]^x + ^I + ly) (9) 

iV/o 4 

2 x,fl<u 
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Figure 2: Phase diagram obtained from the MC simulation. 


where Ni^ = 6 V, Ni^ = 12V and Ni^ = 4V. 

With these dehnitions i?o 0 when j3 ^ oo and > 1 when Ki —> oo. 

On the three-dimensional (/5, ki, K 2 ) parameter space we consider the 
plane (3 = 2.3. On this plane there is a PT line {k\, We expect to learn 
on the properties of this PT on the region 7 ^ 0, 7 ^ 0), where the signals 

are clearer, with regards to applying what we learn to the standard case 

K2 = 0 . 

To monitorize the strength of the phase transition, we measure the exis¬ 
tence of latent heat, and the behavior of the specihc heat. As we shall see, for 
Ki = 0 the transition is first order, with a clearly measnrable latent heat. We 
will see how this transition weakens along the PT line for increasing valnes 
of Ki- 


3 Numerical study 

We have simulated the model in a lattice with periodic bonndary con¬ 
ditions. For the update we have employed a combination of heat-bath and 
over-relaxation algorithms (ten over-relax sweeps followed by a heat-bath 
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Figure 3: Normalized energy distribution for L=12, A = oo, ki= 0.395, /3=2.3 from 
[p^ (dotted line), compared with the distribution we obtain at the same couplings 
(solid line) in L=12 too, when statistics is increased by one order of magnitude. 


sweep). For the simulation we used the RTNN machine, consisting of a net¬ 
work of 32 PentiumPro 200MHz. processor. The total CPU time employed 
has been the equivalent of 3 years of PentiumPro. 

Monte Carlo methods provide information about the thermodynamic 
quantities at a particular value of the couplings. We have used the Spec¬ 
tral Density Method (SDM) [|T^ to extract information on the values of the 
observables in a hnite region around the simulation point. In particular it is 
useful to have a precise location of the coupling where some observables have 
a maximum, as well as an accurate measure of the value of that maximum. 

From the Monte Carlo simulation at some coupling k, we got the his¬ 
togram H (E) which is an approximation to the density of states. Using the 
SDM approximation the probability of hnding the system with an energy E 
at a different coupling k,' can be written as: 

F«/(U) oc (10) 

The region of validity of the SDM approximation is Ak ~ l/(Ucr), being 
a the width of the distribution H{E). Although a gets the maximum values 















in the critical region, the approximation has been very useful, specially for 
tuning the couplings where to measure. 

Concerning the lattice sizes, we have used lattices ranging from L=6 to 
L=24. For the small lattices {L=6, 8 and 12) we have done 4x lO^r iterations, 
being r the largest autocorrelation time for the energy, which ranges from 
r 10 in L=8 to r ~ 35 in L=24. For the largest lattices, L=20 and 24, we 
run up to lO^r MC iterations. 

The statistical errors are computed with the jackknife method. 


4 Results 


We shall make the discussion with the hrst-neighbors link energy, Ei, but as 
far as the critical behavior is concerned, we could carry out the analysis with 
any of the energies. We remark that an appropriate linear combination of 
El, E 2 and Eq could give slightly more accurate results Ei- 

We have considered hxed values of ki (0, 0.02, 0.1, 0.2 and 0.3) and sought 
the K 2 critical for every line, We have also studied the case K 2 = 0 

varying ki which corresponds to the usual SU(2)-Higgs model. 

The SDM has been used to locate the apparent critical point, dehned 
through the specihc heat behavior. From the specihc heat matrix 


Cl’^iL) = 


dKi 


( 11 ) 


we obtain the best signal for C),’^(L), which can be calculated as (we shall 
omit the superscript from now on) 




( 12 ) 


In a hrst order phase transition, behaves, asymptotically, propor¬ 

tional to the volume, L'^. If the FT is second order, the dominant behavior 
for is which diverges too provided that a > 0. At the up¬ 

per critical dimension a = 0, and one has to go further the leading order. 


appearing logarithmic divergences [21 


As a consequence of this divergent behavior, in a hnite lattice Cy{L) shows 
a peak at some value of the coupling which will be taken as apparent critical 
point, K 2 (-^)- 

In Figure ^ we plot the critical line {kI,K 2 )- This line is obtained by 
extrapolation to the thermodynamic limit according to K^ioo) = ^^(T) — 
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Figure 4: Ei distribution on the axis ki = 0 on a L=8 lattice, at {k 2 = 0.15,/9 = 
2.3). The cosine fit is very accurate in spite of the finite (5 value, in this region of 
parameters the pure gauge term couples slightly both sub-lattices. 


AL~^. We point out that this extrapolation is valid only in hrst order phase 
transitions. If the PT is second order the power in L is (—1/z/). 

In continuous PT scale invariance holds, and the thermodynamic magni¬ 
tudes, such as the specihc heat, or susceptibilities do scale. However if the 
transition is hrst order the correlation length remains hnite and hence there 
is not scaling properties, and no critical exponents can be dehned. 

Nevertheless in a hrst order PT we can ask how large is the lattice we 
need in order to observe the asymptotic behavior of C^. In particular, with 
abuse of language one can measure a pseudo ajv exponent to get insight on 
the nature of the PT: the larger is the lattice we need to measure v = 1/d, the 
weaker the PT is. Following this, hrst order PT can be classihed according 
to their degree of weakness. 

However, the so called weak hrst-order PT appear often in literature (see 


^1 and references there in) as PT characterized by a transient behavior 
with a non-measurable latent heat. Let be the correlation length of the 
system at the critical point in the thermodynamic limit. In a hnite lattice 
of size L, the hrst order behavior will be evidenced if L > For lattice 
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Figure 5: Normalized distribution of Ei at (ki = 0 .2, K 2 = 0.10036) in a L=16 
lattice. The cubic fit at the maxima to get AE is superimposed. 


sizes much smaller than the system will behave like in a second-order PT, 
since the correlation length is effectively inhnite. As an example, in Figure 
^ we plot the energy distribution in a L=12 lattice at (/d = 2.3, Ki = 0.395) 
obtained in , compared with the one we obtain at the same couplings and 


in the same L, when the statistics increases by one order of magnitude. We 
observe that the order of the PT can not be discerned at this volume, even 
when the statistics is enough. Termalization effects can also contribute to 
mistake the histogram structure. 

The entire line (^ 1 ,^ 2 ) is hrst order, but the weak character increases 
as ^2 —*■ 0. We will make a quantitative description of the weakening phe¬ 
nomenon by studying the specihc heat, and the latent heat. 

But before going on, we shall make a remark concerning the behavior 
on Ki = 0. As we pointed out, the system is symmetric under the change 
Ki —*■ —Ki- The transformation (^) maps the positive Ki semi-plane with 
energy Ei, onto the negative Ki semi-plane with energy -Ei. The transition 
across this axis is hrst order because the energy is discontinuous. In the limit 
/3 —> cx), and in = 0 , the system decouples in two independent sublattices, 
each one constituted by the hrst neighbors of the other. The hrst neighbors 
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MC iterations 

Figure 6: MC evolution of Ei at (ki = 0.2, for L=8 (lower part), L=12 

(middle) and L=16 (upper part). 


energy for this system is proportional to cos 6, being 9 the angle between 
the symmetry breaking direction of the scalar held in both sub-lattices. In 
Figure ^ we plot the Ei distribution for L =8 at (ki = 0,^2 = 0.15) and 
f3 = 2.3. We see that the agreement with a cosine distribution is quite good 
in spite of the hnite (3 value. 

4.1 Latent Heat 

Along the apparent critical line we have done simulations for different lattice 
sizes and stored the plaquette and links energies to construct the histograms 
for the energy distributions. In a hrst-order phase transition the energy has 
a discontinuity which manifest in the appearance of latent heat, /S.E. This 
quantity is not well dehned in a hnite lattice, so we measure the distance 
between the two maximum of the energy distribution, and extrapolate to 
the thermodynamic limit. The drawback of this approximation is that the 
maxima of the energy distribution are difficult to discern, since this function 
at the apparent critical point is very noisy We have used a cubic spline at 
the maxima in order to get a more reliable estimation (Figure 5). 
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In Figure we show the MC evolution of Ei for L=8, 12 and 16 at 
Ki = 0.2. In L=8 the latent heat is not clearly measurable. We observe in 
the MC evolution how the two-state signal becomes cleaner as the lattice size 
increases (see Figure |^). 

In Figure ^ we plot the distribution of Ei at ki = 0.2 at the apparent 
critical point K^iL) for L=6, 8, 12 and 16. A remarkable stability of AEi 
with the volume is observed. This is a common feature for all values of Ki. 

As we have already pointed out, the transition weakens when increasing 
Ki and larger lattices are needed in order to observe a measurable latent 
heat. We give a quantitative description of this fact in Figure |^, where the 
distribution of Ei for several values of {ki, is displayed. The two- 

state signal is no longer measurable in L=12 at ki = 0.3. The hrst evidences 
of two-states appear in L=20 (see Figure |). but from its energy distribution 
we can only give an approximate value for AEi{L = 20) since the two-peaks 
appear too close to each other. 

In Figure ^ we plot AEi{L) and AE 2 {L) as a function of 1/L^, in order 
to get AEi in the thermodynamic limit with a linear £t. Finally we quote 
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Figure 8: Normalized distributions of Ei for ki = 0.02, 0.1, 0.2 and 0.3 at K 2 ( 12 ). 


these values in Table |l], together with the change in the action (§) between 
the two phases. 

From the energy distributions at fi :2 = 0, see Figure ^ we have no direct 
evidences of the existence of latent heat. However, on the larger lattices one 
can observe non-gaussianities in the energy distributions. Such asymmetries 
could precede the onset of clear two-peak structures in larger lattices, however 
this is just a guess. We conclude that no information concerning the order 
of the FT can be obtained from the energy distributions up to L=24. 

4.2 Specific Heat 

We have done MC simulations at the points predicted by SDM, (ki, ^^(T)), 
in order to measure accurately the peak of Cy{L). 

As an example we show in Figure ^the value of Cy{L) around its maxi¬ 
mum for various lattice sizes at Ki = 0.2 (upper plane) and at Ki = 0.3 (lower 
plane). We observe that the maximum of the specific heat grows slower when 
increasing ki, indicating a weakening in the FT. 

In Figure |T^ we draw C™®'^(L) relative to C™'^^( 6 ) as a function of the 
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Figure 9: Normalized distribution of Ei at ki = 0.3 for L=12,16 and 20 


lattice size. The values have been normalized to C™'^^(8)/C™'^^(6) in order 
to compare distinctly the behaviors for different ki. The slope of the seg¬ 
ment joining the values of in consecutive lattices gives the pseudo ajv 
exponent. We observe that such slope is approximately 4 at Ui = 0.02, 0.1 
and 0.2 for all the volumes we compare. However the transition at ki = 0.3 
evidences much more weakness. We do not have evidences of asymptoticity 
in till L=20, as could be expected from the energy distributions (Fig¬ 
ure H). The slope of the segment joining C™®'^(16) with C™®'^(20) is 3.05(12) 
which is almost the asymptotic value expected for a hrst order FT. 

As expected, only if the two-peak structure is observed in the energy 
distributions and AF^ is stable, the maximum of the specihc heat will grow 
up like the volume, L'^. 

At fi ;2 = 0 we are within the transient region even for L=24. We remark 
that in this case is defined by the element of the specific heat matrix 
(|TT|) since this is the most natural choice at this point, and also is the best 
signal we measure. As we observe in Figure at ^2 = 0 seems to 

tend to a constant value as V—>■ cxo up to L=20. From our previous discussions 
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0 0.01 0.02 0.03 

1/L^ 

Figure 10: Ai^i (upper plane) and AE 2 (lower plane) as a function of l/L"^. The 
two-peak structure is not clearly observed in L=6 at any ki value. The values 
quoted for this lattice size are upper bounds. 


we should conclude that either the correlation length at the transition point 
is much larger than the lattice size up to L=20, or the transition is second 
order with a = 0 in the thermodynamic limit. However, in L = 24 things 
are changing, C™“^(L = 24) starts to run away of this quasi-plateau, and the 
pseudo ajv exponent grows again. As can be observed in Figure |^, at 
= 0.3 the pseudo ajv exponent decreases for the segment L = 12-16, with 
respect to the value in the segment L = 8-12. The lattice L = 20 is enough 
to overcome the transient region, but the behavior is qualitatively the same 
that in K 2 = 0 , though the transition is stronger. 

We believe that this behavior is general for weak first order transitions 
in four dimensions. There exists a transient region in which the correlation 
length is effectively infinite compared with the lattice size, and the system 
behaves like suffering a second order FT with thermal index a 0 in the 
thermodynamic limit. 
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Figure 11: Normalized distribution of Ei at K 2 = 0 , in L=16, 20 and 24 at the 
peak of the specific heat 


4.3 Binder Cumulant 

In order to check the consistency of our results we have also considered the 
behavior of the Binder cumulant 


This quantity behaves differently depending on the order of the PT. If the 
transition is second order the minimum of the cumulant, approach 2/3 
in the thermodynamic limit. However if the transition is hrst order, 
tends a value smaller than 2/3 reflecting the non-gaussianity of the energy 
distribution at the transition point. 

In Figure |T^ we plot for several ki values. 

For those values of Ki in which the PT is distinctly hrst order, the mini¬ 
mum of the Binder cumulant stays safely away from 2/3, as we observe 
extrapolated to L—>• oo is 0.65401(4) at ki = 0.1. However, this value reachs 
0.66637(5) at Ki = 0.3, and 0.66657(8) at K 2 = 0. Again we hnd a tight 
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coupling 

AE'i(cx)) 

AE2(oo) 

A^ 

= 0 

- 

0.0366(8) 

0.0134(12) 

= 0.02 

0.0162(6) 

0.0347(5) 

0.0137(13) 

Ki = 0.1 

0.0162(7) 

0.0345(9) 

0.0094(10) 

Ki = 0.2 

0.0179(7) 

0.0201(8) 

0.0078(12) 

Ki = 0.3 

f«0.006 

^ 0.012 

0.0026 


Table 1: AE(oo) for Ei and E 2 , and variation of the action. 


difference between a very weak first order PT and a continuous one. 

For the sake of discussing quantitatively the order of magnitude of the 
latent heat in the limit K 2 = 0, from the energy distributions we find that 
at K 2 = 0 one can approximately locate one the peaks of the energy at Eq 
0.273. The other should be at certain E;, = — A, being A the latent heat. 

In the thermodynamic limit the energy distributions are two delta func¬ 
tions situated at E^ and E^, then 


2(g* + Ej) 
HEl + 


( 14 ) 


If we use V™" = 0.66657 and Eq in (pA|) the value we got for the latent 
heat is A 0.006 which is of the same order as the one expected from the 
histograms. 


5 Conclusions 

The order of the Confinement-Higgs phase transition in the SU(2)-Higgs 
model with fixed modulus is a highly non trivial issue. We have used an 
extended parameter space, in order to get a global vision on the problem. 
On this extended parameter space we have found a line of hrst order phase 
transitions which get weaker as K 2 —> 0. We have also observed that, because 
of the computer resources needed, it is too ambitious trying to measure two- 
peak energy distributions in the limit ^2 = 0. However, on this point, we 
can get conclusions from the behavior of the specihc heat. 

As we have discussed along the paper, a fake second order PT seems 
to exists for a range of L in very weak hrst order phase transitions. We 
have applied Finite Size Scaling properties along this transient region to 
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Figure 12: at ni = 0.2 (lower plane), at ki = 0.3 (middle) and at K 2 = 0 

(upper plane). The dotted line is the SDM extrapolation. 


compute a pseudo ajv critical index. We want to be extremely careful at 
this point, this computation is completely meaningless when the transition 
is first order, since Scaling does not hold, but it can be used as a technical 
tool to catalogue the PT when there is no direct evidences, as in this case. 
Using the relation a = 2 — i/d, we got v varying in the interval (0.36, 0.41) 
in the range L = 8 ,..., 20. Calculated from L=20 and L=24, v k, 0.35. We 
expect this behavior to be transitory, and when going to larger lattices sizes, 
if the transition is second order, v should reach its mean held value v = 1 / 2 . 
If the transition is hrst order this value should go to 1/d, indicating that 
the specihc heat maximum grows like the volume We believe that this 
is the case, since the pseudo u exponent in L=24, instead of approaching 
1/2, starts to decrease. An example of weak hrst order PT showing a similar 
behavior is described in ||25|| . 

In what concerning the motivation of introducing a second coupling, we 
pointed out that K 2 should not change the order of the PT because do not 
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In (L/6) 

Figure 13: for the various ni values and K 2 = 0. We have normalized the 

values with respect to C™®'^(8)/C™'^^(6). The slope of the segments is indicated 
when smaller than 4. 


change the symmetry properties. This argument is heuristic, but the phase 
diagram we found supports this assertion. As far as the order of the PT is 
concerned, we think that this approach can be useful when dealing with PT 
of questionable order in the sense that it is not clear whether the transition 
is weakly hrst order or higher order. The hope is that it could be applied to 
other more controversial models. 

I thank A. Tarancon and L.A. Fernandez for comments and advice. 
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